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1. Four and five flavors 

As has been known for a long time [1,2] processes relying on gluon splitting into an initial state 
bb pair can be described by generating a universal bottom content in the proton via the DGLAP 
evolution and including it as a partem density. In our case this means that to leading order we 
compute the process P5P5 — > tH with five quark flavors inside the proton [3, 4, 5, 6]. Following 
Figure 1 we can as well limit the quark content of the proton to four light flavors and only generate 
initial state bottom quarks via gluon splitting, so our leading order process is P4P4 — > btH~ [7]. 
The index of the proton indicates that this implies different proton pictures, i.e. different parton 
densities and a s values. 

To all orders in perturbation theory the two approaches are equivalent. The difference between 
them is the ordering in perturbative QCD. The DGLAP equation resums collinear logarithms which 
appear at every fixed order in a s . Since they diverge for massless quarks there is no question that for 
those the logarithms have to be resummed, to obtain finite and stable rate predictions. For bottoms 
the logarithms have the form log(M/m/ ; ) with the process dependent hard scale M ~ ra, + ra#. They 
are regularized and reside in the perturbative regime nit, » Aqcd- 

Because of these logarithms, comparisons at leading order have to be taken with a grain of 
salt. The renormalization and factorization scale dependences in both processes are large and 
the predicted rates correspondingly uncertain. In the five-flavor scheme the factorization scale 
dependence of the total rate is approximately logarithmic and has a one-to-one correspondence 
to the collinear logarithm [1]. Setting the factorization scale to /if = M implicitly assumes the 
existence of large logarithms all the way to log(M/m&). A careful study of the kinematics of gluon 
splitting generating the logarithms shows that bottom densities combined with an incoming gluon 
or bottom never generate as large logarithms, due to the steep drop of the gluon parton densities 
towards large x [2, 4, 8]. Only a consistent (lower) choice of the bottom factorization scale ensures a 
stable perturbative behavior of the five-flavor production rate. For tH production these effects are 
not huge, while they completely change the balance of the different production modes for example 
in the case of charged Higgs pair production [9]. Moreover, given that we have fully understood 
NLO results for the five-flavor [4, 5, 6] and four-flavor [7] descriptions of tH" production we do 
not expect a huge discrepancy between the different rate predictions. 

Expanding the argument to distributions there appear slight differences between the two cal- 
culations. The distributions of the heavy top and Higgs are computed with the same perturbative 
accuracy as the respective total rates. Luckily, we know for example from the production of super- 
symmetric particles [10] that the distributions of heavy particles produced in pairs are remarkably 
stable with respect to higher orders in QCD, so again we do not expect large discrepancies. 




Figure 1: Feynman diagrams for associated top-Higgs production at the LHC. 
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Figure 2: Transverse momentum distribution of the top-Higgs pair, simulated at fixed order (NLO, five 
flavors) and including the parton shower [6]. The assumed charged Higgs mass is 300 GeV. 

2. MC@NLO and jet radiation 

The main difference between the NLO four-flavor and five-flavor calculations appears in the 
bottom jet distributions. On the one hand, the motivation for using bottom partons is that in most 
events we will not see the forward and comparably soft bottom jet, and integrating over its phase 
space produces the logarithmic rate enhancement. On the other hand, modern analysis methods 
are able to make use of such features which are only visible in a fraction of the events, so our 
simulation has to describe them correctly. 

To leading order, in the five-flavor scheme the heavy pair recoils against a strictly longitudinal 
bottom. Only at NLO this recoiling jet acquires a finite pj : b [5]. In the four-flavor scheme the 
recoiling bottom is realistically described already at leading order. At first glance this looks like 
a case for a four-flavor calculation, at least when looking at distributions. However, from the 
moderate size of the collinear logarithm we know that pjj, *C M, i.e. in the bulk of its phase space 
the bottom jet is correctly described by a parton shower and not by a hard matrix element. What we 
should really do is combine either the four-flavor or the five-flavor description with a parton shower. 
This is precisely what MC@NLO does. The impact of the parton shower becomes obvious in the 
px distribution of the top-Higgs system shown in Figure 2. At fixed order it diverges while the 
resummed collinear logarithms in the parton shower correctly describe the low-pr regime. 

Combining the parton shower with the matrix element to NLO ensures that we correctly de- 
scribe all jet radiation while the total rate and the hard distributions are at NLO accuracy. This way 
we do not need an external normalization, as it used to be the case for PYTHIA and HERWIG. 
With a built-in NLO rate prediction we can answer questions of the kind [6]: 

• What are the transverse momentum or rapidity spectra of the radiated b jet? 

• Are these b spectra the same as for light-flavor initial state radiation? 

• Is the b jet from the top decay always harder than the b jet from initial state radiation? 

• How likely is it that we observe or tag this radiated b jet? 

• How often do we see light-flavor or b jets given one b jet from the top decay has been observed? 
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Table 1: Probability (%) to see a b or light jet in a given detector region, after one jet, most likely from the 
top decay, is observed. The top decays leptonically or hadronically while the charged Higgs is stable [6]. 



As an example, we answer the last question in Table 1 [6]. The starting point, motivated by 
the usual analyses, is that one bottom jet from the top decay should always be observed. The 
additional bottom jet arises from initial state gluon splitting. In the five-flavor scheme to leading 
order the probability of seeing such a bottom jet within the fiducial volume of the detector is zero. 
Including the NLO corrections with a reliable pr t, spectrum over the entire range indicates that in 
roughly 18% of the events we can tag the radiated bottom jet given pjj, > 25 GeV and \ r]b\ < 2.5 
for run = 300 GeV. Combined with typical tagging efficiencies this translates into roughly 10% of 
the inclusive tH events with a b tag. Because of the collinear enhancement the radiated bottom 
jet becomes harder for larger Higgs masses, but the generic picture of roughly 10% of all events 
including a tagged forward bottom does not change. 

In addition, we can ask if the appearance of a forward jet can be viewed as a sign for initial- 
state gluon splitting into two bottoms. Table 1 indicates that for a leptonic top decay, i.e. no 
additional hard jets in the event, it is two to three times as likely to observe a light-flavor jet from 
initial state or final state radiation than a bottom jet. When we in addition allow for jets from a 
hadronic top decay almost all events will see an additional jet with pjj > 25 GeV and \r\j\ < 2.5. 
This percentile only drops once pjj reaches the Jacobian peak of the top decay. 

Again, only the combination of a hard NLO matrix element with the parton shower allows us 
to systematically compare such kinematic features of bottom jet radiation, light-flavor jet radiation, 
and top decay jets, as shown in Table. 1 [6]. 

3. Bottom mass and bottom partons 

Bottom parton densities are based on the splitting of an off-shell gluon into a pair of massive 
bottoms. While for light-flavor quarks the splitting threshold is of the order of Aqcd and hence not 
numerically relevant, for bottoms it is in the range of perturbative QCD. This makes it a relevant 
input parameter in the computation of bottom parton densities, which has to be tested [1 1]. To first 
approximation, a shift in the bottom mass changes the logarithmic parton densities by 
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Table 2: Absolute and relative production rates for tH production at NLO, varying the input bottom mass in 
the on-shell scheme. The coupling is fixed by tan j3 = 30, and the renormalization scale is jx = [m t + m#)/2. 



M AT M ( 8m b \ M 8m h 
log >\og — — =log log H ^log . (3.1) 



' Jll}, nib + $ m b m b \ m b ) m b m b 

For well motivated applications of the bottom parton densities (M S> nib) the uncertainly due to the 
bottom mass becomes increasingly irrelevant. For example, a logarithm of the order log 100 ~ 4.6 
reduces the relative error from the bottom mass input by this factor 1 /4.6. 

Because in the matrix element squared of the hard process the bottom mass only enters power 
suppressed we set it to zero in that part of the calculation. This includes the collinear divergence 
which we regularize in dimensional regularization and move into the parton densities. It is dis- 
putable if together with the input bottom mass, i.e. the threshold of the gluon splitting, we should 
also vary the bottom Yukawa coupling. However, for large enough tan/3 the change of the produc- 
tion rate with the bottom Yukawa coupling is trivially given by a oc y*. If we evaluate the bottom 
Yukawa in the MS scheme, changing the on-shell bottom mass by ±0.5 GeV changes yb{M) by 
roughly ±13% and the Higgs cross section by ±25%. 1 

In comparison, varying the bottom mass in the parton densities over 4.25 •••5.25 GeV for 
7 TeV collider energy leads to a ±5% window for a Higgs mass of 200 GeV, as shown in Table 2. 
The change in rate is half of the change in the bottom mass, roughly in line with Eq.(3.1). For 
larger collider energies and for larger Higgs masses the relative change in rate slightly decreases, 
again as expected. This again confirms Eq.(3.1) and its conclusion that the effect of the bottom 
mass through the parton densities is strongly suppressed in comparison with the shift in the bottom 
Yukawa. 



4. Small Higgs masses 

For small Higgs masses a problem arises in the NLO corrections to the process bg — > tH : for 
example the 0(yfaj) subprocess 

gg/qq^tH-b (4.1) 

includes diagrams with on-shell gg — > it production and a subsequent decay i — > H b. They only 
include a two-particle production phase space, a t f x BR f -_^ ff -^, and therefore numerically dominate. 

'Note that in the MC@NLO implementation the starting value m2,(m&) = 4.23 GeV is fixed, corresponding to an 
on-shell value of 4.87 GeV. Varying the input bottom mass will not actually change the bottom Yukawa coupling. 
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Figure 3: Total cross sections for top-Higgs production at the LHC including small Higgs masses. Left: 
PROSPINO implementation [5]; Right: MC@NLO implementation [6], for the correct diagram subtraction 
and the inconsistent diagram removal. 



From a field theory perspective there exists no problem with these on-shell diagrams, aside 
from the usual complications of intermediate on-shell particles. To compute LHC rates and distri- 
butions we can introduce a width Y t or apply more advances schemes in analogy to four fermion 
production at LEP [12]. From Tevatron top analyses we expect that the introduction of a width 
while keeping the spin correlation of the top quark will be completely sufficient. 

The two-fold reason for separating top pairs from top-Higgs production is phenomenological: 
first, the two processes have different kinematic features and therefore require different analyses 
based on different (simulated) event samples. Second, the QCD corrections to the two processes 
should not be expected to be identical. Therefore, we first simulate top pair production including 
QCD effects and normalized to the measured rate. Top-Higgs production then corresponds to an 
event sample which we can add to the top pair sample without any double counting. 

The appropriate treatment of such intermediate states in different 2 — > 2 production channels 
has been developed in PROSPINO [10] and the NLO pp ->■ H t calculation [5]. It was then 
independently applied to single top production in MC@NLO as so-called diagram subtraction [13]. 
The challenge is to define local subtraction terms which remove the on-shell process pp — » tt from 
the real emission contribution Eq.(4.1). Interference between the two processes of course counts 
towards new physics, i.e. top-Higgs production. This is why diagram removal schemes on the 
amplitude level are not useful. 

In terms of the momentum squared Sbh flowing through the intermediate top propagator and 
the corresponding residue |^#| (sbh) we add a subtraction term to the otherwise unmodified Breit- 
Wigner top propagator 

\^\ 2 (sbh) \^£\ 2 (m}) 9n . 

7 2,2 7 2r2 " 7 2ul 2t* ®^~^t) ®{m t -m H ) , (4.2) 

[s bh -mfy + mfTf [s bh -mf) l + mfTf 

including the appropriate constrained phase space [10]. This expression holds independently of 
the value of the regulator T t . In the PROSPINO implementation with quasi-stable supersymmetric 
particles we use the gauge invariant narrow width limit to uniquely define the subtraction term 

r 1 ™ 7 m f\l\ 2 r 2 =7lS ( S bh ~ m ]) ■ ( 4 -3) 
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Table 3: NLO Production rates in pb for pp — > tH production in the five-flavor scheme, using the in- 
put parameters defined by the LHC Higgs cross sections working group. The central renormalization and 
factorization scales are chosen as jio = (mt + mn) /4 [4, 5]. The bottom Yukawa scales with tan/3 = 30. 

In contrast, the MC@NLO implementation includes top decays, so we identify the regulator with 
the physical anomalous top width. 

In Figure 3 we show the different contributions to associated top-Higgs production at the 
LHC [5]. Top pair production including a Breit-Wigner propagator and the appropriate branching 
ratio is labelled (Jew- It dominates for small Higgs masses and quickly becomes numerically irrele- 
vant in the threshold region. This description might be improved by including all off-shell diagrams 
contributing to tH~b production. However, <7 ff-sheii om y sn o ws an effect above threshold, where 
it corresponds to a subset of the NLO top-Higgs diagrams. Instead, we use the complete subtracted 
Oinci,NLO which together with the on-shell top pair production gives the full rate a sum . The resulting 
K factor for the subtracted tH production is flat across the threshold m# = m t . 

As a side remark, the four-flavor prediction for p — > btH~ [7] includes the on-shell top con- 
tribution to leading order, regularized by the physical top width. However, the separation into 
two distinctly different processes would work the same way, for example if one would want to 
normalize the top-pair production rate to data. 

5. Outlook 

We have reviewed the status of theoretical predictions for the production rate of associated 
tH production in a two-Higgs doublet model. Two schemes describe bottom initiated processes 
at hadron colliders, with or without bottoms densities inside the proton. To all orders in perturbative 
QCD they are equivalent, but to leading order they give different results. Both schemes have been 
implemented to NLO, and we find reasonable agreement for the production rates. 

For distributions the difference between these two schemes is not the main issue. The relevant 
collinear bottom and light-flavor jet radiation is be poorly described by either of them and requires a 
combination with a parton shower. MC@NLO describes the kinematics of bottom and light-flavor 
jet radiation reliably over the entire phase space and allows us to study the features of radiated 
bottom jets, radiated light-flavor jets, and jets from the top decay. 

One yet unstudied issue in the five-flavor scheme is the dependence of the rate on the bottom 
mass assumed for the bottom parton densities. The effect of a shift in the input bottom mass is 
suppressed by the same collinear logarithm which justifies the use of bottom parton densities. Its 
numerical impact does not exceed a few per-cent. 
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For light charged Higgs bosons it is crucial that we consistently combine higher order correc- 
tions of the kind pp — > tH b with top pair production followed by an anomalous decay t — > H b. 
There exists a unique prescription how to separate the two processes using a point-wise phase 
space subtraction. We have reviewed the basic features of this PROSPINO subtraction scheme and 
its MC@NLO implementation. 

Finally, we give the NLO cross section predictions for the five-flavor scheme using the input 
parameters defined by the LHC Higgs cross section working group. For two different LHC collider 
energies and a range of renormalization and factorization scales they are shown in Table 3. 

We would like to thank the organizers of this workshop series for the very inspiring and con- 
structive atmosphere. Moreover, we would like to thank Michael Kramer for many constructive 
discussions. 
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